Effective medium approach for stiff polymer networks with flexible cross-links 
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Recent experiments have demonstrated that the nonhnear elasticity of in vitro networks of the 
biopolymer actin is dramatically altered in the presence of a flexible cross-linker such as the abundant 
cytoskeletal protein filamin. The basic principles of such networks remain poorly understood. Here 
we describe an effective medium theory of flexibly cross-linked stiff polymer networks. We argue 
that the response of the cross-links can be fully attributed to entropic stiffening, while softening due 
to domain unfolding can be ignored. The network is modeled as a collection of randomly oriented 
rods connected by flexible cross-links to an elastic continuum. This effective medium is treated in 
a linear elastic limit as well as in a more general framework, in which the medium self-consistently 
represents the nonlinear network behavior. This model predicts that the nonlinear elastic response 
sets in at strains proportional to cross-linker length and inversely proportional to filament length. 
Furthermore, we find that the differential modulus scales linearly with the stress in the stiffening 
regime. These results are in excellent agreement with bulk rheology data. 

PACS numbers: 87.16.Ka, 87.15. La, 82.35.Pq 



I. INTRODUCTION 



The mechanical response and locomotion of living cells 
is mainly controlled by the cellular cytoskeleton. The 
cytoskeleton is a highly composite network of various 
stiff biopolymers, along with various binding proteins 
for force generation, cross-linking and polymer growth 
regulation. Understanding the basic physics that gov- 
erns the mechanical properties of a composite biopoly- 
mer network represents an important biophysical chal- 
lenge that will help elucidate the mechanics of a liv- 
ing cell. In addition to their importance for cell me- 
chanics, cytoskeletal networks have also demonstrated 
novel rheological properties, especially in numerous in 
vitro studies [E S, S i, II H, 0, i, H [H- However, 
there have been few theoretical or experimental stud- 
ies that address the coniposite nature of the cytoskele- 
ton [m, [mil El, [H, Recent experiments on F- 
actin networks with the highly compliant cross-linker fil- 
amin, in particular, have demonstrated several striking 
features: These networks can have a linear modulus as 
low as 1 Pa, which is significantly lower than for actin gels 
with incompliant cross-links, and yet they can withstand 
stresses of 100 Pa or more and can stiffen dramatically by 
up to a factor of 1000 under applied shear [10, 11]. Both 
the linear and nonlinear elastic properties of actin-filamin 
gels appear to be dramatically affected by the flexible 
nature of the cross-links, resulting in novel behavior as 
compared to actin-networks with incompliant cross- links, 
and to synthetic polymer gels. This suggests new network 
design principles that may be extended to novel synthetic 



materials with engineered cross-links (12] . However, the 
basic physics of networks with flexible cross-links remain 
unclear. 

In this article we provide a detailed description of an 
effective medium approach to describe the nonlinear elas- 
tic properties of composite networks consistin g of stiff 
filaments linked by highly flexible cross-links 
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schematic image of the network we aim to model is shown 
in Fig. [T] The network is composed of randomly oriented 
filaments/rods of length L, which are linked together by 
highly flexible cross-linkers. The cross-links consist of 
two binding domains interconnected by a thermally fluc- 
tuating flexible polymer chain of length £o- The compli- 
ance of such a cross-linker is entropic in nature. 

Adopting the WLC model, we can fully characterize 
the cross-linkers with a contour length and a persis- 
tence length [13, 111]- The WLC force-extension curve, 
which is shown in Fig. [5] c) demonstrates the dramatic 
stiffening of the cross-linker as it reaches its full ex- 
tension. Indeed, atomic force microscope (AFM) mea- 
surements show that an actin cross-linker such as fil- 
amin can be accurately described as a wormlike chain 
(WLC) [H, [13]. At large mechanical loads, however, the 
experimental force-extension curve deviates significantly 
from WLC behavior. The polymer chain in cross-linkers 
such as filamin consists of repeated folded protein do- 
mains, which unfold reversibly at sufficiently large me- 
chanical loads. The experiments by Furuike et al. [l^] 
show that after an initial stiffening regime at a force- 
threshold of « 100 pN one of the protein domains un- 
folds reversibly. The accompanied increase in contour 
length results in a strong decrease in the cross-linkers 
stiffness. This softening is immediately followed by WLC 
stiffening as the thermal undulations of the lengthened 
cross-linker are stretched out. This leads to an elastic 
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response that alternates between entropic stiffening and 
softening caused by domain unfolding, resulting in a saw- 
tooth force-extension curve. 

It has been suggested that the unfolding behavior of 
filamin is crucial for the mechanical properties of net- 
works with such cross-linkers [lli ITsl l20j|. Simulations 
of stiff polymer networks, assuming a sawtooth force- 
extension curve for the unfoldable cross-links, reveal that 
such networks exhibit a fragile state in which a signifi- 
cant fraction of cross-linkers is at the threshold of domain 
unfolding [13] . This results in strain softening of the net- 
work under shear, inconsistent with the pronounced stiff- 
enin g respo nse observed experimentally in actin-filamin 
gels We estimate, however, that under typical in 

vitro experimental conditions, domain unfolding in the 
cross-links is highly unlikely. For domain unfolding to 
occur with multiple filamin crosslinks with tensions of 
order 100 pN, the resulting tension in the actin filaments 
is likely to exceed rupture forces of order 300 pN of F- 
actin [2l|. Also, a simple estimate of the macroscopic 
stress corresponding to even a small fraction of filamins 
under 100 pN tensions is larger than the typical limit of 
shear stress before network failure is observed. Therefore, 
we do not expect domain unfolding to occur. Rather, it 
seems likely that cross-link unbinding occurs before suffi- 
ciently large sufhciently large forces are attained for a sig- 
nificant amount of domain unfolding. Detailed estimates 
based on experiments suggest filamin tensions only of or- 
der 1-5 pN at network failure UM. It has also been shown 
in single molecule experiments [22| that filamin unbinds 
from F-actin at forces well below the forces required for 
unfolding, which indicates that cross-linker unfolding is 
highly unlikely to occur in typical network conditions. 
Therefore, we consider only the initial stiffening of the 
cross-links, which we show can account well for the ob- 
served nonlinear elasticity of actin-filamin gels. 

Our model consists of a network of stiff filaments 
connected by flexible cross-linkers. The compliance of 
such a network is expected to be governed by the cross- 
linkers. The stiff filaments provide connectivity to the 
network and constraint the deformation of the cross- 
linkers, thereby setting the length scale of the effective 
unit cell of the network. Consequently, we expect that 
the elasticity of the network will be controlled by the 
filament length L and network connectivity, which is ex- 
pressed in terms of the number of cross-link per filament 
n. Therefore, we describe the network with a model in 
which the basic elastic element consists of a single stiff 
rod and many compliant cross-linkers that are connected 
to a surrounding linear elastic medium. 



II. EFFECTIVE MEDIUM APPROACH 

Networks of semiflexible polymers with point-like in- 
compliant cross-links have been studied extensively [3, [^, 
m, [3i m, HEl • These systems exhibit two distinct elastic 
regimes: One in which the deformation is affine down to 




FIG. 1: Schematic figure of an isotropic stiff polymer network 
with highly compliant cross-linkers. 



the smallest length scales of the network and a regime 
that is characterized by highly non-affine deformations. 
Simulations [2X, 2Sl have shown that the deformation of 
these networks becomes more affine with increasing cross- 
link concentration and polymer length length, which has 
been borne out by experiments [4, 29]. The elastic re- 
sponse of the network can fully be accounted for by the 
stretching modes of the polymers in the affine regime. 
In addition to stretching modes, stiff polymers can also 
store energy in a non-affine bending mode. Indeed, it has 
been shown that in sparser networks, in which there are 
fewer constraint on the constituting polymers, non-affine 
bending modes dominate the elastic response [U, [l^, . 
We will, however, not consider the sparse network limit 
here. 

We expect the soft stretching modes of the cross-linkers 
to govern the elasticity of a dense network of stiff poly- 
mers with highly flexible cross-links. However, the large 
separation in size and stiffness between cross-links and 
filaments does imply a non-uniform deformation field for 
the cross-links at the sub-filament level. On a coarse- 
grained level the network deforms affinely and stretches 
the cross-links as depicted in Fig. [5] b). The network 
surrounding this particular rod is shown here as a grey 
background. The deformation of the cross-links increases 
linearly from in the center towards a maximum value 
at the boundaries of the rod. At small strains the cross- 
links are very soft and follow the deformation of the stiffer 
surrounding medium. However, at a strain ~ ^o/L 
the outer-most cross-links reach their full extension and, 
consequently, stiffen dramatically. This suggest the ex- 
istence of a characteristic strain 7c, for the onset of the 
nonlinear response of the network. 

The macroscopic elasticity of the network results from 
the tensions in all the constituting filaments. The tension 
in a particular filament can be determined by summing 
up the forces exerted by the cross-links on one side of the 
midpoint of the filament. We will employ an effective 
medium approach to calculate these forces as a function 
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of filament orientation and the macroscopic strain. Thus, 
we model the network surrounding one particular rod, as 
a continuum, which effectively represents the elasticity of 
the network, as depicted in Fig. [2] a) and b). We then 
proceed by considering contributions from rods over all 
orientations to calculate the macroscopic response of the 
network. 

The remainder of this article is organized as follows. 
First we study a model in which the effective medium is 
treated as a linear elastic continuum. In this model we 
will describe the cross-links both as linear springs with fi- 
nite extension, and also as WLC cross-links. We analyze 
our model in both a fully 3D network, as well as a simpli- 
fied ID representation, which already captures the essen- 
tial physics of the nonlinear behavior. At large strains, 
when many of the cross-linkers are extended well into 
their nonlinear regimes, it is no longer realistic to model 
the surrounding network as a linear medium. Therefore, 
we extend our linear medium model in a self-consistent 
manner, replacing the embedding medium by a nonlin- 
ear effective medium whose elastic properties are deter- 
mined by those of the constituent rods and linkers. This 
self-consistent model can quantitatively account for the 
nonlinear response found in prior experiments on actin 
filamin networks [Hill- Finally, we show how we can 
compute the tension profiles along the filaments and we 
demonstrate how to use these to express the macroscopic 
stress in terms of the maximum force experienced by a 
single cross-link. 



III. THE LINEAR MEDIUM MODEL 

We first develop a one dimensional representation of 
our model, which will be used in section |V] to construct 
a more realistic three dimensional model. Also we will 
restrict the treatment here to a linear description of the 
effective medium, a constrain that we lift in section ITVl 

Consider a rigid rod of length L connected by n flexible 
cross-links to an elastic medium. We shall refer to such an 
elastic unit as a Hairy Rod (HR). The medium is subject 
to an externally imposed extensional strain e parallel to 
the orientation of the rod. The presence of the HR in 
the medium reduces the deformation of the medium at 
a position x in the rest frame of the rod by an amount 
UEAiix, e) = ex — Uci{x, e), where Uci{x, e) is the extension 
of a cross-linker at a distance x from the center of the 
rod. The magnitude of Uci{x, e) and uem{x, e) are set by 
requiring force balance between the cross-links and the 
medium. 

fci{uci{x, e)) = KemUem{x, e), (1) 

where fd (u) is the force-extension curve of a single cross- 
linker. The tension tq in the center of the rod is found 
by summing up the forces exerted by the stretched cross- 
links on one side of the midpoint of the rod. Assuming 
a high, uniform line density n/L of cross-links along the 
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FIG. 2: (Color online) a) a single filament connected by n flex- 
ible cross-links to the surrounding network, which we model 
as an efi'ective elastic continuum (shown here as a grey back- 
ground) and b) illustrates the proposed nonuniform defor- 
mation of the cross-linkers on a single filament in a sheared 
background medium, c) Force-extension curve of a Hookean 
Finite Extendable (HFE) cross-linker (dashed blue curve) and 
of a WLC cross-linker (solid black curve) . 

rod, we can write the sum as an integral 

'^o(e) = y / dx' fci{uci{x' ,e)). (2) 
^ Jo 

where Uci{x',e) is found by solving Eqn. ([1]). The full 
tension profile r(e, x) is found by replacing the lower limit 
of the integration by x 

T{e,x) = —J dx' fci{uci{x' ,e)) (3) 

A. Hookean finite extendable cross-linkers 

We can solve Eqns. ([T]) and ^ to compute the the mid- 
point tension in a rod, as soon as a force-extension curve 
for the cross- links is specified. In the absence of unfolding 
or unbinding, we can describe the force-extension behav- 
ior of a flexible cross-linker such as filamin with the WLC 
model, as depicted with the black solid line in Fig.[2]c). It 
is instructive to simplify the WLC force-extension curve 
by assuming a Hookean response with a spring constant 
kci up to an extension £q, which is the molecular weight 
of the cross-linker. The spring constant kd = is 
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found from the WLC model for small extensions in the 
limit £p <§; £q [l^, where fc^T is the thermal energy. Be- 
yond an extension Iq, the cross-linker becomes infinitely 
stiff. The force-extension curve of these Hookean Finite 
Extendable (HFE) cross-links is shown as a blue dashed 
curve in Fig. [2]c). The finite extensibility of the cross- 



links implies a critical strain €c 



L/2 



at which the cross- 



linkers at the boundaries of the rod reach full extension. 
For strains e < Ec 



To(e) = Y / dx' 
^ Jo 



kclKEM ^ 

kci + Kem 



(4) 



Thus, the midpoint tension depends linearly on strain for 
e < tc- For larger strains, the expression for the midpoint 
tension in a hairy rod in Eq. ^ reads 



ro(e) 



, kclKEM I 

dx — ex 



+ 



dx' 



kci + Kem 

kclKEM 

kci + Kem 



(5) 



£o + KEMiex' - io) 



The expression has separated into two integrals, clearly 
representing a sum over the cross-links with an extension 
< £o and a sum over the cross-links that have already 
reached full extension. We also note that beyond ec the 
midpoint tension depends nonlinearly on strain. Using 
Eq. ([5]) we compute the ID modulus Gid — tq/c, as 
shown in Fig. [3l Below the critical strain, the response 
is dominated by the linear elasticity of the cross-links 
GiD ~ \nkciL. The cross-links at the edge of the rod 
become rigid at a strain threshold Ec = 2£qIL. As the 
strain is further increased, the outer cross-links stiffen 
consecutively, resulting in a sharp increase of G\d- At 
large strains, Gid asymptotically approaches a second 
linear regime ^ ^uKemL. 



B. Worm Like Chain Cross-Linkers 

We now consider flexible cross-linkers described by the 
more realistic WLC force-extension curve, as depicted by 
the solid line in Fig. [2]c). The force-extension relation is 
well described by the interpolation formula Ui 



fci{u) 



kuT 



( 




(6) 



where fc^T is the thermal energy. Using Eqs. ([T]) and ^ 
we can calculate the ID modulus Gid for cross-linkers 
with this force-extension curve. The result of this cal- 
culation is shown in Fig. |3) The force-extension curve 
of the WLC cross-linker is linear up to extensions very 
close to £q, upon which a pronounced stiffening occurs, as 
shown in Fig.[5]c). We can exploit this, together with the 
property that for a dense network the medium is much 
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FIG. 3: (Color online) a) The modulus Gizj = ro/e for the ID 
representation of the linear medium model with HFE cross- 
links with Kem ~ Wkd (blue dashed dotted curve) and 
Kem = lOOfec; (red dotted curve). We also show Gid for 
the model with WLC cross- links with Kem = Wkd (blue 
dashed curve) and Kem = lOOfcci (red solid curve). The in- 
set shows the ratio of the average tension f and the midpoint 
tension tq. 



stiffer than the flexible cross-linkers Kem ^ kd to write 
an approximate expression for the tension in a hairy rod 
in a closed form analogous to Eq. ([5]). 



'ro(e) = T / dx' [ 
^ Jo Jo 



du- 



kcl{u)KEM 

kci{u) + Kem 

kcl{u)KEM 



(7) 



du 

kcl[u) + KEM 



+ Kem {ex' - £o) 



where kci{u) is the differential stiffness dfd/du of the 
WLC cross-linker. This equation states that an HR unit 
deforms essentially affine up to the critical strain. Be- 
yond ec, those cross-links that have reached full exten- 
sion are no longer compliant and start to pull back on 
the surrounding medium. The approximate calculation 
of dTQ/dj using Eq. ([7]) is shown together with the ex- 
act calculation performed with Eq. ([2]) in Fig. |4l This 
graph demonstrates that the approximation captures the 
essential behavior of the exact curve, and results only in 
a minor quantitative difference in the cross-over regime. 
Therefore, we will continue constructing our model using 
this approximation. 

The ID modulus calculated with Eq. ([7]) is shown for 
the WLC cross-links together with the results of the HFE 
cross-links in Fig. [3l Although the main behavior is very 
similar to that of the HFE cross-linker model, the use 
of the more realistic WLC force-extension curve has in- 
troduced a considerable smoothing of the cross-over. The 
nonlinear behavior in the WLC force-extension curve ini- 
tiates slowly well before full extension, resulting in a more 
gradual onset of nonlinear behavior of the HR with WLC 
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FIG. 4: (Color online) The ID modulus Gid = r/e of the 
rod with WLC cross-linkers as a function of the extensional 
strain e imposed on the medium parallel to the orientation 
of the rod. The solid red curve shows the exact calculation 
using Eqs. ^ and ([2]) and the dashed blue curve shows the 
approximate calculation using Eq. (O. 



cross-linkers. Remarkably, the characteristic strain ec for 
the nonlinear behavior is proportional to (q/L, indepen- 
dent of the exact nonlinear response of the linkers. 

For a calculation of network mechanics the average ten- 
sion f in a filament is more relevant than the midpoint 
tension [soj. f is found by averaging the tension pro- 
file given by Eq. ([3]) along the backbone of the filament. 
The ratio t/tq is shown in the inset of Fig. [3] We find 
that over a broad range of strains f = 3/2to. During 
the cross-over regime the ratio exhibits a peak with an 
amplitude that depends on the exact ratio of Kem and 
kci- 



IV. SELF-CONSISTENT MEDIUM MODEL 

The linear treatment of the effective medium breaks 
down at large strains. The network, consisting of a collec- 
tion of many HR's, will exhibit nonlinear response when 
the cross-linkers start to get extended into their nonlin- 
ear regime. Thus, it is no longer realistic to assume that 
the effective medium, which should reflect the network 
elasticity, remains linear. In this section we extend our 
model by requiring that the elasticity of the background 
medium self- consistently represents the nonlinear elastic- 
ity of the constituent HR's. The elasticity of the medium 
should therefore depend on the density of filaments and 
on the elasticity of an HR averaged over all orientations. 
Thus, we require the stiffness per cross-link of the effec- 
tive medium Kem to be determined by the stiffness of 
an HR 
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a dro 
nL de 



(8) 



The proportionality constant a depends on the detailed 
structure of the network. In section IIV Al we derive an 
expression for a in the continuum elastic limit. The mid- 
point tension in a rod can be written down analogous to 
Eq. (0) 



n 

roie) = - 



dx' : 



de' 



kci{x'e') 



a. dr I X f 







(9) 



where kci{u) is the derivative of the force-extension re- 
lation of the cross-linker. Note that we have applied 
the same approximation as we did in Eq. ([7]). However, 
we expect this approximation to hold even better here, 
since the medium stiffens strongly as well as the cross- 
links. Eq. ([9]) can be simplified to the following differen- 
tial equation for tq (e) 



' de 
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de2 

4 fe,,(eL/2) + ^±21 

a drp 
4 de 



(10) 



1^ ^ - L/2 



We find the following behavior of the model with WLC 
cross-linkers: Below the characteristic strain for nonlin- 
ear response ec = 2£o/^j the tension in a rod depends 
approximately linearly on strain. This linearity will be 
reflected in the self-consistent effective medium, and con- 
sequently, the model shows behavior similar to the lin- 
ear medium model up to the critical strain. By solving 
Eq. (|10p we find the midpoint tension tq in a rod as a 
function of extensional strain e. Beyond the critical strain 
the tension depends highly nonlinearly on strain, with a 
derivative that increases as 



dTQ 
de 



(11) 



Note that unlike in the linear medium model, where the 
derivative asymptotes to a final value set by Kem, here 
dro/dj increases indefinitely. For the HFE cross- linkers 
we find similar behavior, although in that case the cross- 
over between the linear regime and the asymptotic stiff- 
ening regime is more abrupt. 



A. Continuum elastic limit 

Here we derive an expression for a in the continuum 
elastic limit. Note that this will only be a good approx- 
imation for a dense, isotropic network. The modulus of 
the medium Gem can be expressed in terms of the stiff- 
ness ^ of a HR by averaging over rod orientations [l^l 



r - 1 ^ 

^network _ P , ; 
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(12) 



where p is the length of filament per unit volume, p can 
also be expressed in terms of the mesh-size p = 1/^^. In 
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the linear medium treatment in section Hill we foimd that 
f = |ro. Thus, the network modulus reads 



G 



network 



2 1 dro 
456 W 



(13) 



We proceed by relating Kem to Gnetwork , which enables 
us to find an expression for a. Consider a rigid rod of 
diameter a and length L, which we use as a microrheo- 
logical probe in an effective elastic medium with a shear 
modulus Gem- If the rod is displaced along its axis, 
it will induce a medium deformation 5i that leads to a 
restoring force acting along its backbone. The restoring 
force per unit length is given by 2TrGEM/ ^og{L/a) x Si. 
Here we ignore the log term, which is of order 2tt. Thus, 
the stiffness of the medium per cross-link Kem is related 
to Gem by 



Kem — —Gem- 
n 



(14) 



By requiring Gem ~ Gnetwork we find a from Eqns. (fTB]) 
and 
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(15) 



Note that for a dense network a ^ 1. 



V. 3D NETWORK CALCULATION 

In this section we describe in detail how the macro- 
scopic mechanical properties of a uniformly deforming 
network can be inferred from single filament proper- 
ties. This procedure has been used to describe the vis- 
coelastic [111 and nonlinear elastic properties 0, H, Q of 
semiflexible polymer networks with point-like rigid cross- 
links, although a detailed derivation of this theory is still 
lacking. The main assumption of this calculation is a 
uniform, or affine deformation of the network. The va- 
lidity of the affine treatment of cross-linked semiflexible 
polymer networks has been subject to much debate. In- 
terestingly, 2D simulations in the zero temperature limit 
have found that the deformation can be both affine and 
non-affine depen ding on the density of the network and 
filament rigidity [271 [28j. Here we derive the affine the- 
ory for the case of a filamentous network with point-like 
rigid cross-links. Then we show how this framework can 
be used together with the effective medium approach to 
describe the mechanics of stiff polymer networks with 
fiexible cross-links. 

Consider a segment of a filament between two cross- 
links with an initial orientation h. When subjected to a 
deformation described by the Cauchy deformation ten- 
sor Aij , this filament segment experiences an extensional 
strain directed along its backbone 



This extensional strain leads either to compression or ex- 
tension in the polymer segment depending on its orien- 
tation, and thus results in a tension T(|An| — 1). The 
contribution of this tension to the macroscopic stress de- 
pends also on the orientation of the polymer segment. By 
integrating over contributions of the tension over all ori- 
entations accordingly, we can compute the macroscopic 
stress tensor aij. We calculate the contribution of the 
tension in a polymer segment with an initial orientation 
•h as follows. The deformation A^ transforms the orien- 
tation of the segment into n'j = Ajkni~/\A.n\. Thus, the 
length density of polymers with an orientation n that 
cross the j-plane is given by g^Ajfenfe, where the factor 
detA accounts for the volume change associated with the 
deformation. For the network calculations in this article 
we consider only simple shear, which conserves volume 
(detA =1). The tension in the i-direction in a filament 
with an initial orientation n, as it reorients under strain, 
is r(|An| — ljAiini/\Ah\. Thus, the (symmetric) stress 
tensor reads 



P 



detA 



T{\An\ 



\An\ 



(17) 



The angular brackets indicate an average over the initial 
orientation of the polymer chains. 

One remarkable feature follows directly from Eq. (|17p . 
A nonlinear force extension curve for the filaments is not 
strictly required for a nonlinear network response [3]| . 
To demonstrate this we express the extensional strain of 
a filament explicitly in terms of the strain tensor 7^; 



Vl + 2ukinkni - 1. 



(18) 



Anl 



1. 



(16) 



Thus the extensional strain of a filament depends non- 
linearly on the macroscopic strain of the network. Ad- 
ditionally, the reorientation of the filaments under strain 
leads to an increasingly more anisotropic filament distri- 
bution. Remarkably, these geometric effects result in a 
stiffening of the shear modulus under shear strains of or- 
der 1, even in the case of Hookean filaments. At large 
strains all filaments are effectively oriented in the strain 
direction, which limits the amount of stiffening to a fac- 
tor of 4 (2D networks) and 5 (3D networks) over the 
linear modulus at strains of order 10. Thus the stiffening 
due to this effect occurs only at large strains and is lim- 
ited to a factor 5. Therefore we expect this mechanism 
to have a marginal contribution to the more dramatic 
stiffening that is observed in biopolymer gels at strains 
< 1 0, la| • We would like to stress that the geomet- 
ric stiffening discussed above has a different nature than 
the geometric stiffening discussed by [s^. In their 
case, the stiffening is attributed to a cross-over between 
an elastic response dominated by soft bending modes in 
the zero strain limit and a stiffer stretching mode dom- 
inated regime at finite strains. In the affine calculation 
described here, only stretching modes are considered. 

By limiting ourself to a small strain limit, we can ex- 
clude the geometric stiffening effects discussed above. 



7 



This is instructive, since it allows us to study network 
stiffening due to filament properties alone, and it is a very 
good approximation for most networks since the nonlin- 
ear response typically sets in at strains < 1. For a volume 
conserving deformation (detA = 1) in the the small strain 
limit the stress tensor in Eq. (fT7|) reduces to [s^l 



where S£„ 



1«2 



i^/ip is the total stored length due to 



(19) 



In this limit the geometric stiffening mechanism discussed 
above is absent. Next we show explicitly how to calculate 
the shear stress a^^^ in the z-plane for a network, which 
is sheared in the a;-direction. A filament segment with 
an orientation given by the usual spherical coordinates 9 
and if undergoes an extensional strain 

e = a/I + 27 cos{ip) sm{9) cos{9) + 7^ cos^ (9) - 1 
« 7Cos(v?)sin(6l)cos(6'), (20) 

where we have used a small strain approximation in the 
second line. The tension in this segment contributes to 
the x2-component of the stress tensor through a geomet- 
ric multiplication factor cos((/)) sin(6') cos(6'), where the 
first two terms are due to a projection of the forces in 
the a:-direction and the second term is due to a projec- 
tion of the orientation of the filament into the orientation 
of the z-plane. The stress in the xz-direction is thus given 
by 

(Jrr. = — / / d9dipsm{9){ (21) 



_P_ 

47r 



T [7 cos{(p) sm{9) cos{9)] cos{ip) sin{9) cos{9) } . 

Since we limit ourselves to the small strain limit, we do 
not account for a redistribution of the filament orienta- 
tions by the shear transformation in this equation. 



A. semifexible polymer networks with rigid 
point-like cross-links 

In this section we show how the affine framework can 
be used to compute the elastic response of a network with 
inextensible semiflexible polymers connected by point- 
like rigid cross-links. 

Consider a segment of an inextensible semiflexible 
polymer of length £c between two rigid cross-links in the 
network. Thermal energy induces undulations in the fila- 
ment, which can be stretched out by an applied tension. 
By adopting the WLC model in the semiflexible limit 
4 > 



-pi 



the force-extension relation of this segment has 



been shown to be given implicitly by [22 



(22) 



7r24 ^ n2(n2-f (/,)' 

where is the tension t normalized by the buckling force 

2 

threshold k-^. This relationship can be inverted to ob- 
tain the tension as a function of the extension St. 

r = {.5tl5l^^^) , (23) 



equilibrium fluctuations. This is also the maximum ex- 
tension, which can be found from Eq. (|22p as ^ cx). 
For small extensions 5£ this reduces to 



90 



Si. 



(24) 



This result can be inserted into Eq. to find the 
linear modulus of the network 



Go = 6/9 



(25) 



For a network in either two or three dimensions, the max- 
imally strained filaments under shear are oriented at a 
45 degree angle with respect to the shear plane, meaning 
that the maximum shear strain is 



7n 



14 

3 £7} 



(26) 



Using the small strain approximation(as in Eq. p9|) ). we 
can calculate the nonlinear network response 



a 

(Jr. 



d9d(p sm{9){ 



(27) 



47r7o JO 

cos{(p) sin{9) cos(6')] cos{ip) sin{9) cos{9)] 

where we define the critical stress to be tic = Ptt- We 
have also defined 7 = 7/7C, where the critical strain for 
the network is 



_ 14 
64' 



(28) 



This equation demonstrates that the nonlinear response 
of a network of inextensible semiflexible polymers with 
rigid cross-links is universal _4] for small strains. Note, 
however, that this would not hold if we would use the full 
nonlinear theory from Eq. p7p , valid for arbitrarily large 
strains. Thus, geometric stiffening effects may lead to 
small departures from universality. Alternatively, univer- 
sality may break down as a result of enthalpic stretching 
of the polymer backbone [5]. 

In this section we have assumed that at zero strain 
all filament segments are at their equilibrium zero-force 
length. However, cross-linking of thermally fluctuat- 
ing polymers will result in cross-linking distances both 
smaller and greater than their equilibrium length. This 
effect, which is ignored in our discussion here, leads to 
internal stresses build into the network during the gela- 
tion [|. 

The universal nonlinear elastic response for a semiflex- 
ible polymer network with rigid cross-links is shown in 
Figs. [HI and [5| The divergence of the differential modulus 
beyond the critical strain is of the form 



(l-7max)^ ' 

depicted in Fig. [5| This results into a powerlaw stiffening 
regime of the form K ~ a^^"^, as shown in the inset of 
Fig. [S| This prediction is consistent with experiments on 
actin gels with the rigid cross- linker scruin [^]. 
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B. stiff polymer networks with highly flexible 
cross-links 

For a network with flexible cross-links we do not con- 
sider the tension in filament segments, but rather the 
average tension f in the whole filament. By using the 
effective medium approach we can compute the average 
tension in a filament as a function of the orientation of 
the rod and the macroscopic shear strain 7. Contribu- 
tions to the stress from the average tension in the rods 
are integrated over all orientations according to Eq. (PT|) . 
In our description we thus assume affine deformation of 
the network on length scales > L. Note, however, that 
we do not assume that the cross-links deform affinely. 

We find both from the linear medium model and the 
self-consistent model for a network with highly flexi- 
ble cross-links that the linear modulus is approximately 
given by 

Go « -pnkciL. (29) 
o 

The appearance of the filament length L in this equation 
is remarkable, and is due to the non-uniform deformation 
profile of the cross-links, which enhances the forces ap- 
plied by the cross-links further from the midpoint of the 
filament. The onset of nonlinear elastic response occurs 
at a critical strain 

7c = (30) 

The full nonlinear response as predicted by our model 
is shown in Figs. [5] and [51 The results of the linear 
medium model with WLC cross-links, as shown with a 
green dotted line, are qualitatively similar to the results 
of the ID model(see Fig. [3]). For the self-consistent model 
we find that beyond 7c the differential modulus increases 
as a powerlaw, as shown in Fig. [5l Interestingly, we find 
only a small quantitative difference between the model 
with HFE and WLC cross-links. 

The differential modulus K = da/d-f is plotted as a 
function of stress in Fig. [S] The stress is normalized by 
the critical stress CTc, which we define here as 

o-c = Gojc = -pnkciio- (31) 

We find a sharp increase in stiffness beyond the critical 
stress, which quickly asymptotes to a powerlaw regime, 
where the exponent is given by 1 — 1/(^(L/^)^ — 1). In- 
terestingly, this exponent does not depend on the exact 
form of the nonlinear response of the cross-linkers. This 
exponent emerges as a consequence of the finite extend- 
ability of the cross-links and the non-uniform deforma- 
tion profile along the backbone of the filament. Remark- 
ably, the powerlaw exponent is not universal. However, 
in the dense limit we consider in our model, the devia- 
tion to an exponent of 1 is <C 1 and depends only weakly 
on the ratio L/S.. As an example, we consider a typical 
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FIG. 5: (Color online) The differential modulus K — daf/dy 
normalized by the linear modulus Go as a function of strain 
normalized by the critical strain 7c. The universal curve for 
a semiflexble polymer network with rigid cross-links is shown 
as a black dashed curve. We also show the results of the 
self-consistent model with WLC cross-links (red solid curve) 
and simple cross-links (blue dash-dotted curve), the linear 
medium model with WLC cross-links with Keai = WOkd 
(green dotted curve) 

in vitro network for which ^ = 0.3^m and the average 
filament length is L = 15/im. For this case we find an 
exponent of 0.98. The asymptotic powerlaw regime with 
an exponent « 1, as predicted by our model is consistent 
with recent experimental data on actin networks cross- 
linked by filamin [ll|, [l^ . 

The inset of Fig. [H] shows the rigid linker model to- 
gether with the self-consistent model for a network with 
flexible cross- links. In this case the stress is normalized 
by a stress cto, which marks the knee of the curve. 



VI. TENSION PROFILES AND SINGLE 
CROSS-LINKER FORCE ESTIMATE 

Recently, there has been much debate on the mechan- 
ical response of actin binding proteins such as filamin. 
Specifically, it is discussed whether the cross-links stiffen, 
unfold or unbind under tension in both physiological or 
in vitro conditions. This issue has major implications 
for the dynamical and mechanical properties of the cy- 
toskeleton. The discussion has been partially resolved 
recently by single molecule [22] and bulk rheology [l6j 
experiments on the actin-filamin system. These experi- 
ments indicate that cross-links unbind at forces well be- 
low the force required for domain unfolding. It is crucial 
for the bulk rheology experiment, to be able to infer the 
forces experienced by a single cross-linker from the mea- 
sured mechanical stress. In this section we show that by 
using the shape of the tension profile, we can relate a 
macroscopic quantity such as the stress to the maximum 
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in the linear response regime Eq. (fT7|) simplifies to 
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FIG. 6: (Color online) The differential modulus K = da/ /d'y 
normalized by the linear modulus Go as a function of stress 
normalized by the critical stress CTc for the self-consistent 
model with WLC cross-links (red solid curve), HFE cross- 
hnks (blue deish-dotted curve) and the linear medium model 
with WLC cross-links with Kem = lOOfcci (green dotted 
curve). The inset shows the rigid linker model together with 
the self-consistent model for a network with flexible cross- 
links. In this case the stress is normalized by a stress (Jq, 
which marks the knee of the curve. 



force experienced by a single cross-linker in the network. 

The tension along a single filament is not uniform in 
networks of stiff finite length filaments and incompliant 
cross-links [13, HJ]. It was found in simulations that 
in the afRne regime the tension profile is flat close to 
the midpoint and the tension decreases exponentially 
towards the boundaries of the filament. In the non- 
affine regime a different tension profile has been re- 
ported, in which the tension decreases linearly towards 
the ends [s^l- In the case of a flexibly cross-linked net- 
work of stiff polymers we also expect a non-uniform ten- 
sion profile, although in this case the underlying physics 
is different. The deformation of a cross- linker at a dis- 
tance X from the midpoint of the rod is Ud x"f and, con- 
sequentially, cross-links further away from the midpoint 
exert larger forces on the rod, resulting in a non-uniform 
tension profile. 

We can calculate the tension profile for a given rod 
using Eq. In the limit of highly flexible cross- linkers, 
the tension profile in the linear elastic regime is given by 



L kci + Kem 2 




(32) 



The tension profiles as computed with the self-consistent 
model with WLC cross-links are shown for various strains 
in Fig. [71 For low strains we find a parabolic profile, 
which flattens out towards the edges for larger strains. 

We now proceed to estimate the force experienced by 
a single cross-linker. For an affinely deforming network 



(33) 



Filaments at a 45° angle with respect to the stress plane 
bear the largest tension fmax and experience a strain 
along their backbone of 7/2. Assuming linear response 
we find fmaxij) = T'i'y)/'^- In the case of a parabolic ten- 
sion profile, the average tension f in a filament is related 
to the largest force /o experienced by a cross-linker at the 
boundary of the rod by t = ^n/g. Thus, we can express 
the macroscopic stress in terms of the maximum forces 
experienced by cross-linkers on the filaments under the 
greatest load 



^ — , ^ f max • 

45 



(34) 



For the derivation of this equation we have assumed to 
be in the linear response regime. In the nonlinear regime 
we expect the expression to still hold approximately, al- 
though the prefactors will change. 

Kasza et al. [1^1 found that the failure stress of the net- 
work Umax is proportional to the number of cross-links 
per filament n in actin networks with the flexible cross- 
linker filamin. This suggests that filamin failure, rather 
than rupture of single actin filaments is the cause for 
network breakage. In contrast, for actin networks with 
the rigid cross-linker scruin, which binds more strongly 
to actin than filamin, rupture of actin was found to be 
the mechanism for network failure On the basis of 
our model and the experimental data from Ref. [l6j we 
estimate filamin failure forces of order 1 — 5 pN, far be- 
low the unfolding force 100 pN. This suggests that net- 
work failure is due to filamin unbinding. This is con- 
sistent with recent single molecule experiments, which 
show that filamin unbinding is favored over unfolding of 
the Ig-domains for low loading rates [2^ . 



VII. IMPLICATIONS AND DISCUSSION 

We have studied the nonlinear elasticity of stiff poly- 
mer network with highly flexible cross- links. We find that 
the mechanics of such a network is controlled by network 
connectivity expressed in the number of cross-links per 
filament n. This was found earlier in experiments on 
actin-filamin gels , providing strong experimental ev- 
idence for cross-link dominated mechanics in these net- 
works. Within this picture, stiffening occurs at a strain 
where the cross-links are stretched towards their full ex- 
tension. As a result, we expect 7c to be proportional to 
the molecular weight of the cross-linker £o- This predic- 
tion is consistent with the results of Wagner et al. [l3] , 
where cross-link length was varied, while keeping the av- 
erage filament length fixed. Interestingly, they observed 
larger values of 7c than expected either from our model 
or based on Refs [O, [H, [li . 
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FIG. 7: (Color online) The reduced tension profile along the 
rod, normalized by the midpoint tension tq. This profile is 
calculated with the self-consistent model with WLC cross- 
hnks. 

In addition, we find here that the filament length L 
plays an important role in the nonlinear elasticity of these 
networks. In particular, the onset of nonlinear response 
7c ^ ^o/L depends crucially on filament length. This has 
been confirmed by recent experiments on actin-filamin 
gels, showing an approximate inverse dependence of the 
7c on actin filament length [s^. The sensitivity of net- 
work response to filament length, both in experiments 
and in our model, appears to be one of the hallmarks of 
actin-filamin networks. On the one hand, this may ex- 
plain the apparent difference between the critical strains 
reported in Refs. [H [H, [3. On the other hand, this 
also suggests that it may be more even more important in 
such flexibly cross-linked networks to directly control and 
measure the filament length distribution than for other 
in vitro actin studies ^Z^. Our model does not account 
for filament length polydispersity. A distribution in fila- 
ment length is expected to smooth somewhat the sharp 
stiffening transition predicted by our model. 

The dependence of the critical strain for networks with 
flexible cross- links observed in experiments and predicted 
by our model is in striking contrast with the behavior 
found for rigidly cross- linked networks. In the latter case 
theory predicts 7c ~ ^p/^c (see Eq. (^51) ). which is con- 
sistent with experimental observations The insensi- 
tivity of the nonlinear elasticity of dense networks cross- 
linked with rigid linkers to filament length would suggest 
that network mechanics cannot be effectively controlled 
by actin polymerization regulation. We have shown here 
that the filament length plays a crucial role for networks 
with flexible cross-links, which are abundant in the cel- 
lular cytoskeleton. Thus regulating actin length by bind- 
ing/capping proteins such as gelsolin may enable the cell 
not only to sensitively tune the linear elastic modulus, 
but also the onset of the nonlinear response of its cy- 
toskeleton. 



In the nonlinear regime we expect the differential mod- 
ulus to increase linearly with stress for a dense flexibly 
cross-linked network. This behavior is a direct conse- 
quence of the non-uniform deformation profile along a 
filament and the flnite cxtendability of the cross-links, 
although it is independent of the exact shape of the 
force-extension behavior of the cross-links. The powerlaw 
stiffening K with y « 1 is consistent with recent 

experiments on actin-filamin gels [ll|, [lB| ■ This stiffen- 
ing behavior is very different from the nonlinear response 
observed for actin gels with rigid cross-links for which a 
powerlaw exponent of 3/2 is observed [4j], consistent with 
theory for an affine response governed by the stretching 
out of thermal fluctuations of the actin filaments. Inter- 
estingly, in vivo experiments show that cells also exhibit 
power-law stiffening with an exponent of 1 [svt . 

In this article we have examined a limit in which the 
stiffness of the cross-links is small compared to the stiff- 
ness of an F-actin segment between adjacent cross-links. 
For a large flexible cross-linker such as fllamin this is 
clearly a good approximation in the linear regime. How- 
ever, as the cross-links stiffen strongly they could, in prin- 
ciple, become as stiff as the actin segment. This would 
have a dramatic consequence for the nonlinear response 
of the network. To investigate this we have calculated 
the differential stiffness df/ du as a function of force / for 
a fllamin cross-linker and an actin segment with a length 
0.5 to 2 ^m, spanning the range of typical distances 
between cross-links in dense and sparse networks respec- 
tively. This result is shown in Fig. [51 We flnd the dif- 
ferential stiffness of a fllamin cross-link is always smaller 
than for an F-actin segment, even at large forces in the 
nonlinear regime. This justifies our approach, in which 
we have ignored the compliance of the actin, for a broad 
range of experimentally accessible polymer/cross-linking 
densities. However, at sufficiently high filamin concentra- 
tions, it may be possible that individual network nodes 
involve multiple cross- linkers, in which case the actin fil- 
ament compliance may also become relevant. Thus the 
affect of the compliance of F-actin remains an interesting 
topic for further research. 

We also use our model to study these networks on a 
more microscopic level, such as the non-uniform tension 
profiles along the filament backbone. These profiles can 
be used to establish a relation between the macroscopic 
stress and the largest force experienced by a single cross- 
linker in the network. This allows us to estimate the 
forces experienced by filamin cross-links under typical in 
vitro and in vivo conditions. We find that the load on 
these cross-links is not sufficiently high to lead to signif- 
icant domain unfolding of the filamin Ig-domains, even 
at stresses large enough to rupture the network. Indeed 
both rheology experiments on actin filamin gels and sin- 
gle molecule experiments indicate that unbinding occurs 
well before domain unfolding. 

In other large flexible cross- links such as spectrin [s^ . 
domain unfolding occurs at lower, more relevant forces. 
In this case the domain unfolding could have a dramatic 
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FIG. 8: (Color online) The differential stiffness df/du as a 
function of force / for a filamin cross-linker (solid line) and for 
several F-actin polymer segment lengths. 



affect on the nonlinear viscoelasticity of such networks. 
In previous work, DiDonna and Levine have simulated 2D 
cross-linked networks, where they have assumed a saw- 
tooth force-extension curve for the cross-linkers to mimic 
domain unfolding .1^ . They report a fragile state with 
shear softening when an appreciable number of cross- 
linkers are at the threshold of domain unfolding. Our 
model is based on the stiffening of the cross-linkers, which 
occurs at forces far below those required for domain un- 



folding. This leads to strain stiffening at a point where 
only a fraction of cross-linkers are at their threshold for 
nonlinear response. Thus in both our model and that of 
Ref. [Ts^ the network responds strongly to small strain 
changes, though in an opposite manner: stiffening in the 
present case vs softening in Ref. [l^ . 

In related work, Dalhaimer, Discher, and Lubensky 
show that isotropic networks linked by large compliant 
cross-linkers exhibit a shear induced ordering transition 
to a nematic phase [T3] . It would be interesting to investi- 
gate the affect of the nonlinear behavior of the cross- links 
on this transition. In the present calculation we have 
assumed an isotropic network. An ordering transition, 
which results in a strong alignment of filaments will dra- 
matically affect the nonlinear elasticity of the network. 

In this article we have studied networks of stiff poly- 
mers linked by highly flexible cross-links. Both exper- 
iments [nl, [3] and our model [31 indicate that these 
networks have novel nonlinear rheological properties. We 
find that the network mechanics is highly tunable. By 
varying filament length, cross-linker length and network 
connectivity we can sensitively regulate the linear and 
nonlinear elasticity over orders of magnitude. These 
unique properties can be exploited in the design of novel 
synthetic materials. 
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